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Abstract. We have conducted an optical study of the 
environments of seven radio-loud quasars at redshifts 1 < 
z < 1.6. In this paper we describe deep B and R band 
images obtained for fields of ~6x6 arcmin around these 
quasars with 3cr limiting magnitudes of B ~ 26 and R ~ 
25.4. We found an statistically significant excess of faint 
(B > 22.5 and R > 22.0) galaxies on scales of r <170" 
and r <35" around the quasars. The number of excess 
galaxies, their magnitudes and the angular extension of 
the excesses are compatible with clusters of galaxies at 
the redshift of the QSOs. We also found that the quasars 
of our sample are in general not located at the peak of 
the density distribution, lying in a projected distance of 
about ~70" from it. This result agrees with the scenario 
of a link between overdensities of galaxies around radio 
sources and the origin of their radio emission. 
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1. Introduction 

Evidence has been accumulating over the last decade 
that there is a link between quasar activity and the 
environment of the host galaxy (e.g., Stockton 1982; Yee 
& Green 1984, 1987; Gehren et al. 1984; Hutchings et al 
1984; Hintzen 1984; Yee 1987; lEllingson et al. 1991 



Hintzen et al. 1991 



Yee & Ellingson 1993 



Fisher et al. 1996;Yamada et al. 1997; Hall et al. 1998 



Hall & Green 1998). In particular, the galaxy-galaxy 



interaction appears to play a substantial role as the 
triggering or fuelling mechanism of the nuclear activity. 
The large fraction of host galaxies (HG) of radio-loud 
quasars undergoing tidal interactions/merging processes 
(e.g. Disney et al. 1995; Hutchings & Neff 1997), the 
high absolute magnitudes of these galaxies (e.g. Carballo 
et al. 1998, and references therein), and the velocity 
distribution of galaxies around the QSOs (Heckman et al. 



1984; Ellingson et al. 1991 ) lend support to this scenario 
( |Yee fc Ellingson 1993[ ) 
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Furthermore, Yee & Green (1987), Ellingson et al. 
(1991) and Yee & Ellingson (1993) presented strong ev- 
idence of a link between radio emission and clustering of 
galaxies around QSOs. They found that there are signif- 
icant differences between the environments of radio-loud 
and radio-quiet quasars. While radio-loud quasars are as- 
sociated with clusters of galaxies of Abell richness class 0- 
2, radio-quiet quasars inhabit environments of Abell rich- 
ness class poorer than 0. 

These results are based on searches of overdensities 
around low-redshift quasars (z <0.7). At higher redshifts 
the number of studies are scarce, or they are focused in 
one or a few objects (Hintzen et al. 1991; Boyle & Couch 
1993, Hutchings et al. 1993, Yamada et al. 1997, Hall 
& Green 1998). It seems that the dichotomy of radio- 
loud/radio-quiet environment extrapolates to higher red- 
shifts. Hintzen et al. (1991) reported on the observation 
of 16 fields around radio-loud QSOs with 0.9 < z < 1.5. 
Their sample was selected as radio-loud quasars from the 
Veron-Cetty & Veron (1984) catalog. They found a signifi- 
cant excess of galaxies within a radius of 15 arcsec around 
the quasars. However, Boyle & Couch (1993) found no 
excess of galaxies around their sample of 27 radio-quiet 
quasars, at the same redshift. 

Hall et al. (1998) and Hall & Green (1998) have re- 
cently reported on an excess of faint galaxies (K >19) 
around a sample of 31 radio- loud quasars, covering a red- 
shift range z—1-2. Although the sample was observed at 
the optical and NIR wavelengths, the excess was found 
in their AT-band images. The excess occurs on two spatial 
scales: a peak at 8 <35" from the quasars and a smooth 
component down to 9 ^100", which is the largest scale 
sampled by their A"-band images. The magnitudes and 
colors of the excess galaxies are consistent with a popu- 
lation of early- type galaxies at the quasar redshift. This 
excess is stronger for the more radio-powerful and steep- 
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spectrum radio-quasars, and it is consistent with a cluster 
of galaxies with Abell richness class ~ 1.5 ± 1.5. 

It is still an open question whether the most important 
factor affecting the evolution of radio-loud QSOs in clus- 
ters is the absolute density of the environment, the density 
within ^0.5Mpc, or even the conditions in the host galaxy 
and any galaxies interacting with it (Hall & Green 1998). 
It is important to determine if clusters around quasars 
at high-z are really embedded in large-scale overdensities, 
or if they are young, less-peaked, smoothed clusters that 
extend to large scales (^2-3Mpc). Paraphrasing Hall & 
Green (1998) "...larger areas (4'x4') around our sample 
need to be imaged to confirm the large-scale (8 ^100") 
galaxy excess around z=l-2 radio-loud quasars and to de- 
termine the true angular extent of this excess...". 

In this paper we report on the results obtained from 
deep CCD imaging survey of the fields around 7 radio- 
quasars selected from the B3-VLA quasar sample (Vigotti 
et al. 1997), in the redshift range 1 < z < 1.6. The quasars 
are steep-spectrum radio-sources and they have been se- 
lected as the brightest ones of this sample at this z range 
((z) ~ 1.4), using the available optical information when 
the selection was made (Vigotti et al. 1997). New optical 
observations show that it is not true any longer (Carballo 
et al. 1999), but they are still on the range of the bright- 
est ones, having (Mb) ~ —26.6. The size of the images 
(6'x6') will allow us to sample the density of galaxies to 
scales down to ^190". In Sect. 2 we describe the observa- 
tions and data analysis. In Sect. 3 we analyse the possible 
projected clusters around the QSOs, describing some of 
their properties. In Sect. 4 we compare with previous re- 
sults, and discuss about possible interpretations. The con- 
clusions are presented in Sect. 5. Throughout this paper 
we adopt Hq=50 km s _1 Mpc _1 and qo=0.5 . 



2. Observation and data processing 

Observations were carried out during the night of 1997 
March 13 at the Cassegrain focus of the 2.2m Telescope 
at Calar Alto (Spain). B and R standard Johnson's filters 
were used and the pixel scale was 0.533 arcsec-pixel -1 . 
The field of view corresponds to ~6' x 6' in all the images. 
The seeing was ~1.6" , and variations along the night were 
lower than 0.2" . 

An standard observing procedure was used. CCD bias 
were measured taking s. exposure time images. Sky flat- 
fields for each filter were obtained during twilight. Dome- 
flats were also obtained in order to check the sky flat-fields 
stability; however, sky flat-fields have been found to work 
properly and to produce better flat-fielding. In order to 
reduce the cosmic ray events, and obtain a proper cleaning 
of them, the following observing procedure was applied: 
two frames of 900 seconds exposure time at the R band 
and three frames (two of 900 seconds and another one of 
600 seconds) at the B band were obtained for each QSO. 



Table 1. QSOs sample 



B3 Name 


z 


M460 
a 408 


-P40S 


B 


R 


0740+380C 


1.063 


-1.29 


7.88 


18.36+0.09 


17.65+0.03 


0926+388 


1.630 


-1.09 


1.84 


19.68+0.09 


19.22+0.03 


1123+395 


1.470 


-0.82 


1.05 


18.54+0.09 


17.96+0.03 


1148+387 


1.303 


-0.94 


4.08 


17.36+0.09 


16.74+0.03 


1206+439B 


1.400 


-0.85 


14.90 


18.41+0.09 


17.39+0.03 


1315+396 


1.560 


-0.49 


3.87 


18.70+0.09 


18.28+0.03 


1435+383 


1.600 


-0.85 


1.91 


18.23+0.09 


17.83+0.03 



in 10 is W Hz- 



The individual frames were reduced using standard 
tasks in IRAF 1 package. First, the corresponding bias was 
subtracted for each image. A flat-field image was obtained 
for each filter using the average of the sky flats with higher 
signal-to- noise and without saturated pixels. The images 
were corrected using these flat-fields, and the error pro- 
duced by this correction was lower than 0.8%. Cosmic ray 
events in each frame were detected and replaced by the 
average of the neighbours using the COSMICRAYS task 
on IRAF packages. 

The final images, with 1800 and 2400 seconds exposure 
time for the R and B band respectively , were obtained 
by coadding the reduced individual frames for each object 
and filter. 

2.1. Photometric calibration 

Flux calibration for each night and filter was carried out 
by a standard procedure using 28 Landolt faint photomet- 
ric standard stars (Landolt 1992, and references therein), 
that were observed for each filter along the night. The 
extinction coefficient and zero-point were obtained by a 
least-square fitting between the instrumental magnitude 
(—2.5 • log 10 (N c /t), where N c is the number of counts), 
and the air-mass (X). The night was photometric, and 
the rms of these least-square fits were 0.086 and 0.027 for 
the B and R band respectively. 

The QSO magnitudes were measured on the images us- 
ing circular apertures centred at the emission peak. The 
aperture radius was fixed for the seven QSOs to 6.4", and 
it was set to the value at which the intensity of the ob- 
ject reached the background for the QSO with the larger 
seeing image. Table 1 lists a summary of these results, in- 
cluding some useful radio properties of the quasars and 
their redshifts. Carballo et al. (1999) report on the opti- 
cal photometry of 73 B3-VLA quasars. They observed the 
seven objects of our sample, in both the B and R band, 
and quote magnitudes consistent with those reported here, 
within an error of ^0.1 mag. 



1 IRAF is distributed by the NOAO, which is operated by 
AURA, Inc., under contract to the NSF 
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2.2. Data processing and number counts 

We searched for galaxies in the frames using the SEx- 
tractor package (Bertin & Arnouts 199£). This program 
detects and deblends objects using an isophotal thresh- 
old above the local sky and produces a catalogue of their 
properties. This catalogue includes information of posi- 
tion, shape, profile, magnitude and type, for each de- 
tected object. Objects are classified using a neural net- 
work, which takes as input parameters 8 isophotal areas, 
the peak intensity and the FWHM (full width at half max- 
imum) of each object. The seeing, defined as the FWHM 
of the stars, is another input parameter. This network de- 
termines a "star-like-index" , which is, by definition, for 
the galaxies and 1 for the stars. In practice, a value of 
~ 0.7 works fine to separate stars and galaxies. 

We selected a detection threshold of 5 connected pix- 
els (~1.4 arcsec 2 ) above 2cr s k y level per pixel, which guar- 
antees that all the detected objects have signal-to-noise 
larger than 4.5. This detection criterium is similar to those 
reported in previous works. For instance, Hall et al. (1998) 
used a catalog of objects detected with a minimum area 
of 1.9 arcsec 2 and a detection threshold of 5a. We note 
here that for an area of 1.9 arcsec 2 our detection thresh- 
old would be 5.2er. We applied SExtractor to the 7 i?-band 
and 7 5-band images using this detection criterium. Af- 
ter cleaning sources close to the frame borders or near to 
large foreground galaxies, we obtained the final catalogs 
for each field and filter. We found that the mean magni- 
tudes of the faintest detected objects in the different fields 
arc 5=25.58±0.35 and i?=24.80±0.28. These are the 4.5a 
limiting magnitude of our catalogs. 

The 3a limiting magnitude for 5 connected pixels of the 
images can be estimated by measuring the standard devia- 
tion on the sky. This yields mean values of 5j 3 i ^ 1 =26.0±0.4 
and i? 1 3 i ^ 1 =25.4±0.3. This would imply a 4.5cr limiting 
magnitude of B 1 4^ <T =25.6±0.4 and i?f i ^ CT =24.9±0.3, which 
are consistent with the magnitudes of the faintest detected 
objects in the catalogs. 

The 3cx limiting magnitude of the i?-band images is 
similar to values previously reported in recent works. Hall 
& Green (1998) and Hall et al. (1998) reported a 3a lim- 
iting magnitude of r ^25.5 mag, for a minimum detection 
area of 1.9arcsec 2 . Taking into account the relation be- 
tween the Thuan-Gunn r-band and the Johnson i?-band 
magnitudes (r = i? + 0.43 + 0.15(B- V) ; Kent 1985), and 
assuming a range of B — V colours between and 2 (for 
typical elliptical galaxies at < z < 2), the Hall & Green 
3ct limiting magnitude would be R ~24. 8-25.1. These au- 
thors study overdensities of galaxies around radio-quasars 
at 1 < z < 2, whereas our maximum redshift is z ^1.6. 
Therefore we conclude that our images are adequate to 
search for galaxies around QSOs at these redshifts. 

In order to test the validity of our method we have 
made a set of simulations. Seven simulated images have 
been built with the same magnitude-to-counts relation, 



background and noise than the i?-band images of our sam- 
ple. The simulated images were built using the ARTDATA 
tasks in the IRAF package. In each image 1200 galaxies, 
uniformly spatially distributed, and with a number-counts 
distribution following a power-law with a slope of 0.4, 
ranging from i?=18.0 to 25.5, were simualted. We added 
750 stars to each image, with magnitudes in the range 14 
to 25 mag, and similar spatial and flux distribution than 
the galaxies. We have artificially increased the star/galaxy 
number ratio with respect to that observed in our images 
in order to test the galaxy detection and classification at 
faint flux levels. 

The whole procedure was applied over these simulated 
images. The results have been compared with the input 
parameters of the objects. We found that 85% of the 
sources were correctly classified down to R <23.5. The 
percentage of objects (galaxies+stars) misclassified grows 
at fainter magnitudes, being around 30% at R >23.5. 
The fraction of galaxies classified as stars is roughly sim- 
ilar to the fraction of stars classified as galaxies. The in- 
put number-counts distribution of galaxies was recovered 
down to R <23.8 (objects detected up to ~13er). A similar 
result is expected for the objects down to B <24.5 mag. 

2.2.1. Number counts 

A proper determination of the field number-counts and its 
standard deviation is crucial to detect excesses of galax- 
ies and to quantify their significances. Some authors take 
the statistical significance using the poissionian standard 
deviation. But, as have been quoted (Yee & Green 1987, 
and references therein; Benitez et al. 1995), galaxies tend 
to cluster, and the true standard deviation depends not 
only on the number of objects (as in the poissonian case) 
but also on the two-point correlation function. Therefore, 
the standard deviation a(r, N) depends on the number 
of galaxies N and on the angular scale r. This is usu- 
ally expressed as a(r, N) — j(r) ■ vN . The factor j(r) 
takes account of the angular scale, and since there is no 
galaxy-galaxy correlations at any scale, it must tend to 1 
at some large scale (Yee & Green 1987). The limit 7 = 1 
corresponds to a poissonian distribution. Ellingson et al. 
(1991, and references therein) set 7(7") to 1.3 for scales of 
~60 arcsec, whereas Benitez et al. (1995) found a value of 
1.5 for scales as large as ~120 arcsec and Hintzen et al. 
(1991) assumed a poissonian distribution for scales of ~15 
arcsec. We have prefered to make an internal measurement 
of the factor j(r). 

The number of field galaxies in the QSOs images was 
measured directly and then compared with similar results 
available in the literature. Two different, and complemen- 
tary, measurements have been made. First, the number of 
galaxies in the outer part of the images, in an anulus from 
ri n =170" to r out = 185" around the QSO, was measured in 
each image. This radius corresponds to ~ 2.5 Mpc at the 
mean redshift of the QSOs and guarantees that there are 
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no contaminations from possible clusters members. The 
mean galaxy-counts obtained for each QSO field and fil- 
ter defines a density of field galaxies. The measurement 
was repeated for different ranges of magnitude, and we 
obtained a number-counts distribution. This method has 
the advantage that possible contaminations from the puta- 
tive cluster members at these scales have been minimized. 
On the other hand it has the disadvantage that there is 
only one measurement for each field, and therefore, the 
standard deviation is not well determined. The measured 
factor 7 with this method would correspond to 7(170"). 

The second method consists in measuring the number 
of galaxies in a grid of 25 circular, not overlapping, ar- 
eas of 35 arcsec radius (1.1 arcmin 2 ). This angular scale 
corresponds to ~0.5 Mpc at the mean z of the quasars. 
To minimize contamination from possible clusters, only 
areas at distances larger than 70 arcsec from the QSOs 
were considered. This method has the adventage that now 
we have 16 measurements of the number of field galaxies 
for each image, resulting in a total of 112 measurements 
for each band. This provides us with a well determined 
7(35") factor. The major caveat is that it is not possible to 
guarantee that the areas are not contaminated by cluster 
members. Recalling that 7(7") 7^ 1 is due to galaxy-galaxy 
correlations, any clustering would produce an increase of 
this factor. Therefore, any contamination by cluster mem- 
bers tends to increase the measured value of 7(35") and 
to reduce the significance of the possible excess. 

Table 2 lists the number of measured field galaxies 
(we will refer them as expected galaxies throughout this 
paper, N cxp ) within an area of ~1 arcmin 2 (a circular 
area of r < 35") at different magnitude ranges, using the 
two above refereed methods. We also list the computed 
-/Vexp taken all detected galaxies into account (i? a ii and 
i? a ii in Table 2). We have re-normalized the number of 
expected galaxies and its standard deviation to the same 
area assumming that the distribution of field galaxies is 
uniform and that the deviation follows the above refereed 
law (cr=7(r) y^N eKp ) . Table 2 also lists the corresponding 
estimates of the 7 factor at both scales. 

Apparently, there is no significant differences between 
both measurements, at any band and magnitude range. As 
expected, the first method gives a lower value of the total 
expected galaxies (since there is little or no contamina- 
tion from the cluster), while the second method provides 
a better value of the standard deviation. The mean value 
of 7(7-) is 1.79±1.05, with the first method, and 1.08±0.11, 
with the second one. Both values can not be simply com- 
pared, since they have been obtained at different angular 
scales, but the larger uncertainty of the first one is clear. If 
7 is 1.08 at r ~35", at larger angular scales its value must 
be between 1.08 and 1 (the poissonian case). Therefore we 
will use the number counts derived from the first method 
and a standard deviation of 1.08-^/iV CX p for scales larger 
than or equal to 35" . This value of the 7 factor at these 
scales is slightly different from that quoted by Ellingson et 



Table 2. 


Number 


of expected 


galaxies 


within 


1.07arcmin 2 


(r <35") 










N cxp a 




N cxp b 




Mag. range 


a/y/N exp a 


a/ViVcxp " 


B all 


5.88±2.29 


0.95 


6.03±3.27 


1.33±0.26 


B < 23.5 


2.61±3.57 


2.21 


2.70±1.62 


0.99±0.20 


B < 22.5 


1.15±1.57 


1.47 


1.01±1.10 


1.10±0.19 


B < 21.5 


0.28±0.50 


0.95 


0.35±0.59 


1.00±0.13 


B < 20.5 


0.25±0.42 


0.84 


0.17±0.40 


1.00±0.19 


B < 19.5 


0.23±0.41 


0.84 


0.07±0.25 


0.95±0.27 


R all 


8.57±3.78 


1.30 


9.09±3.32 


1.10±0.17 


R<23 


5.08±5.63 


2.50 


5.68±2.33 


0.98±0.26 


R<22 


2.52±6.33 


3.99 


2.54±1.80 


1.13±0.16 


R<21 


1.19±3.63 


3.33 


1.01±1.12 


1.11±0.21 


R<20 


0.64±1.65 


2.06 


0.42±0.75 


1.16±0.22 


R<19 


0.37±0.61 


1.01 


0.22±0.51 


1.09±0.19 



a Measured in an anulus of n n =170" and r ou t~ 185" around 

the QSOs 
b Measured in 16 circular areas of r <35" 



al. (1991) who estimated 7 ~ 1.3. We will be comparing 
our results with those obtained using this value. 

As we mention above, this method guarantees lit- 
tle contamination from possible cluster members, but it 
cannot prevent contamination from large-scale structures. 
Postman et al. (1998) have noticed that such structures 
could expand to distances as large as ~6 Mpc (in our 
adopted cosmology). It is not possible to avoid this prob- 
lem with present data, since it needs a measurement over 
several large-scale structure correlation lengths to exti- 
mate the 'true' number-counts. As in the case of 'cluster 
members', any large-scale structure contamination will in- 
crease the 'real' expected number counts. Therefore, any 
overdensity found would be more likely a lower limit, to 
the extent that the mean density derived from our method 
is slightly biased upward by the outer extent of the corre- 
lated large-scale structure. 

Fig. 1 shows the galaxy number-magnitude relation or 
number-counts distribution for the R and B bands. The 
slopes of the distributions are ~0. 37-0. 36 (within the mag- 
nitude range 19.0< R <23.0) and -0.49-0.40 (within the 
magnitude range 19. 5< B <23.5) respectively. The 100% 
completeness magnitudes, denned as the magnitude were 
the number counts drops from a simple power-law, are 
R ~23 and B ~23.7. These completeness magnitudes cor- 
respond to ~30er detections, i.e. photometric errors near 
to the calibration errors (^0.03 mag). Both the expected 
number of galaxies and the slopes obtained are compat- 
ible with the results of other authors down to this com- 
pleteness magnitudes (Metcalfe et al. 1991, Ellingson et al. 
1991, |Boyle fc Couch 1993J , Benftez 



1991. 
et al 



Hintzen et al. 

1995), with variations smaller than 



20%. These 

small offsets in the number-counts quoted by different au- 
thors are mainly due to the different quality of the images 
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Fig. 1. The R and B-band galaxy number-counts derived 
from the QSO images. 



(e.g., different seeing and photometric conditions), small 
differences in the filter transmissions and the use of dif- 
ferent methods to detect and classify the objects (e.g., 
Hintzen et al. 1991, Kron 1980). 

3. Projected clustering around the QSOs 

Since it is not known whether galaxies in the fields of the 
QSOs arc physically associated with them or not, it is not 
clear which are the optimal magnitude range and angular 
scale to search for possible excesses. Therefore, we have 
preferred to use all the available information, measuring 
the number of galaxies at two different scales around the 
QSOs, r\ <170" and r^ <35", and at different magnitude 
ranges. The number of galaxies found around each quasar 
for each filter band has been added and later compared 
with the expected number of galaxies. The number of ex- 
pected galaxies when adding seven fields was obtained us- 
ing the values listed in Table 2, and were corrected for the 
corresponding areas and number of fields: 

^e 7 /; eWs (n) = 7-(n/35") 2 -iV cxp 

and the corresponding standard deviation is 



7 fields i A _ 



(n) = 7(n) 



A 7 , 



7 fields 



<.\X|> 



(n) 



where r, is either n or r2- 

We have splitted the sample in bins of 1 magnitude 
down to the above defined "100% completeness magni- 
tude" and an additional bin for galaxies fainter than this 



magnitude. It is important to note here that all the galax- 
ies, even galaxies fainter than the 100% completeness mag- 
nitude (R > 23 mag, B > 23.7 mag), were detected up to 
the 4.5a level. Moreover, the number of expected galaxies 
was measured in the same images using the same detec- 
tion/classification procedure as the "excess" galaxies, and, 
therefore, they are affected by the same incompleteness ef- 
fects. We define the relative excess ( 

significance of the excess ( 



N-N c: 



N -N B: , 

Ncxp 

) as na 1 . 



as r\ 



and the 



Table 3 lists the number of detected galaxies at both 
scales, the number of expected galaxies and the signifi- 
cance of the excess for each band. The significance was 
computed using 7 =1.08, and also 7 =1.3 by compari- 
son. These values are listed for any range of magnitude, 
indicated by all, and for the above defined bins of magni- 
tude. We found that there is a significant excess of galaxies 
around the QSOs for both angular scales and filter bands. 
The significance of the excess ranges from 3.56<r for the 
B band and the largest angular scale, to 4.87cx for the B 
band and the smallest angular scale. The excesses are still 
significant if we use 7 = 1.3. The overdensity is detected 
both in the B and R bands at a 3 — 4cr level, and, there- 
fore, the significance of the detection is larger than if it 
were detected only in one of the bands. A quadratic prop- 
agation of the significance of the excesses in both bands 
shows that the detection is significant up to 4.60<r (B+R 
all in Table 3). 

The excess presents a tendency with the magnitude 
range, being significant only for the faintest magnitudes 
(up to — 3<r for R > 22 and B > 22.5), for both bands and 
angular scales. These results indicate that most probably 
the excess is not due to intervening galaxies, since there 
is no significant excess down to R < 21 or B < 22.5 
(see Benitez et al. 1995). The excess galaxies could be 
physically associated with the QSOs: at the mean redshift 
of the QSOs (z ~ 1.4), galaxies with R -23 and B -23.5 
would have Mr. — -20.9 and Mb ~-18.6. A fc+e correction 
typical for an E galaxy has been applied, obtained using 
Bruzual (1983) c-model and the GISSEL code (Bruzual & 
Chariot 1993), assuming a Zf or =5, 1 Gyr of initial burst 
and a Salpeter mass function (Salpeter 1955), with masses 
from 0.1 to 125 Mq. These absolute magnitudes are not 
unusual, and therefore, the galaxies are compatible with 
being at the quasar redshifts (Bromley et al. 1998; Muriel 
et al. 1998). Moreover, recent results (Hall & Green 1998) 
show that galaxies with r >23 mag could be associated 
with quasars at z —1.5. Taking into account the conversion 
between R and r magnitudes, their galaxies should have 
R >22.3 mag, being compatible with our result. 

The spatial extension of the excesses is also compatible 
with this interpretation. Fig. 2 shows the radial distribu- 
tion of the relative excess (77) with respect to each quasar 
(in arcsec), for both the B and R band. The mean radial 
distribution using the seven fields is also shown. The over- 
densities present a peak down to 8 ~35"-40" in five of the 
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Table 3. Number of galaxies within r <170" and r <35" from all the fields 



Mag. range 


N t ot 


N cxp 


nai.3 


nai.os 


N t ot 


N cxp 


ncri.z 


WCTi.08 




1080 




r <170" 




75 




r <35" 




B all 


960.82 


2.96 


3.56 


41.18 


4.05 


4.87 


R all 


1569 


1398.60 


3.51 


4.23 


97 


59.99 


3.67 


3.83 


B+R all 






4.60 


5.53 






5.49 


6.20 


B >23.5 


588 


529.9 


1.94 


2.34 


41 


23.89 


2.69 


3.24 


22.5 < B <23.5 


304 


241.92 


3.07 


3.70 


20 


9.80 


2.51 


3.02 


21.5 < B <22.5 


123 


143.64 


-1.32 


-1.59 


7 


5.57 


0.46 


0.55 


20.5 < B <21.5 


38 


31.92 


0.82 


0.99 


4 


1.37 


1.47 


2.09 


R >23.0 


604 


559.44 


1.45 


1.74 


40 


26.25 


2.06 


2.49 


22.0 <i?<23.0 


524 


423.36 


3.76 


4.53 


31 


17.01 


2.61 


3.15 


21.0 <i?<22.0 


253 


219.24 


1.75 


2.11 


17 


8.82 


2.12 


2.55 


20.0 <_R<21.0 


110 


90.72 


1.56 


1.87 


4 


3.71 


0.12 


0.15 


19.0 < R <20.0 


47 


45.36 


0.19 


0.23 


2 


1.75 


0.15 


0.18 



seven fields (at least in one of the two bands). The peak is 
also visible in the mean distribution of the relative excess. 
The excess drops smoothly from this peak, and remains 
flat until 6 ~T60"-170". If the overdensity is due to clus- 
ters of galaxies at the redshift of the QSOs (z ~ 1.4), the 
peak would extend to ~350 kpc and the cut-off would be 
at about ^1500 kpc. These scales are compatible with the 
physical scales of known clusters at lower redshifts (e.g., 
Dressier 1980, Ellingson et al. 1991). 

The overdensity is produced by ~20-40 galaxies in each 
field (for both angular scales), which corresponds to low 
density clusters (Abell richness class between 0-1, Abell 
1958 ). Summarizing, the magnitudes, angular sizes and 
number of galaxies of the excess indicate that most prob- 
ably we are detecting clusters around the quasars. 

3.1. Dependence of excesses with QSO properties 

The significance of the excess is at least ~4.6<r adding 
the information from both filters and all the fields. This 
indicates that the mean significance of the excess, quasar 
to quasar, must be ~1.7ct. It is interesting to test if the 
field-to-field contribution to the total excess is uniform, or 
it depends on any property of the quasars. 

Table 4 lists the number of detected galaxies at r < 35" 
from the quasars in each quasar field, together with the 
significance of the excess (or defect), for each filter band. 
There are objects where the excess is higher than 4.5cr, but 
others present unsignificant deficits. Moreover, although 
mainly all the quasar fields present a central excess, there 
are large differences in the spatial distribution of the ex- 
cess from quasar to quasar (see Fig. 2). This could indicate 
that there is a dependence of the excess with some prop- 
erty of the quasars. 

In order to investigate this possibility we have studied 
possible correlations between the significance of the excess 
at r < 35" and different properties of the radio quasars. 
No correlation was found between the significance of the 



Table 4. Number of galaxies in each field 



B3 name 


iV B 


nai.os 


JVr 


"■0-1.08 








r <35" 




0740+380C 


5 


-0.34 


19 


3.29 


1148+387 


13 


2.72 


13 


1.40 


1206+395 


18 


4.63 


21 


3.93 


1123+387 


11 


1.95 


12 


1.09 


1315+439B 


11 


1.95 


12 


1.09 


1435+396 


8 


0.81 


7 


-0.50 


0926+383 


9 


1.19 


13 


1.40 


Expected 


5.9 




8.6 





excess, nai.os (of the B and R band and at any angu- 
lar scale), and redshift, absolute magnitude, B — R colour 
or radio spectral index of the QSOs. Fig. 3 shows the ten- 
dency found between this significance and the radio power 
of the quasars at 408MHz. The correlation coefficients be- 
tween both parameters are r xy =0.8 for the R band counts, 
r xy =0.4 for the B band, and r xy =0.7 for the mean of both 
(mean without statistic meaning). This tendency could 
be thought to be due to the large contribution of the B3 
1206+439B field, but, removing it, there is still a ten- 
dency for the R band counts (r xy =0.7). Since our sample 
is short to be statistically significant this tendency has to 
be confirmed over larger samples of radio-quasars, but it 
were expected if there is a connection between radio emis- 



sion and clustering ( Yee & Ellingson 1993, and references 
therein). 

4. Comparison with other results at the same z 
range 



Although there 
tering around 
0.9 



are numerous studies about clus- 
QSOs at redshifts lower than 



(|Yee fc Green 1987 



Ellingson et al. 1991 



Yee fc Ellingson 1993[ [Fisher et al. 1996] , and refer 



p 
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Fig. 2. Radial distribution of the relative excess of galaxies (n) with respect to the QSOs (in arcsec) for each quasar 
field and averaged over all the fields (indicated as mean). The thick and thin solid lines correspond to the objects 
detected in the i?-band and i?-band images, respectively. 



ences therein), the number of studies decreases at higher 
redshifts (Hintzen et al. 1991, Boyle & Couch 1993, 



Hall fc Green 199S ) or they are focused in one or a 
few number of fields (Hutchings et al. 1993). At low 
redshift appears to be clear that radio-loud quasars 
inhabit groups or clusters of galaxies with an Abell 
richness class ranging between and 2, and radio-quiet 
sources inhabit less dense groups of galaxies (e.g., 
Ellingson et al. 1991). These results appear to reproduce 
at high redshift: Hintzen et al. (1991) found a significant 
excess of galaxies around 16 radio-loud QSOs with 
0.9< z <1.5 (32 galaxies with R <23 within 15 arcsec 
when the expectation was 19.3), while Boyle & Couch 



(1993) found no significant excesses around 27 radio-quiet 
QSOs in the same range of redshifts. This dichotomy 
may imply a connection between clustering and radio 
emission ( Yee fc Ellingson 1993 ), which would explain 
the tendency between relative excess and radio-power 
found in our sample. Recently, Hall et al. (1998) and Hall 
& Green (1998) reported their results on optical and NIR 
observations of a sample of 31 2=1-2 radio- loud quasars. 
They found a significant excess of faint (K >19) galaxies 
in these fields, which are compatible with being at the 
quasars redshifts. 

We have found overdensities of galaxies around our 
radio quasars similar to those reported by Hintzen et al. 
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Fig. 2. Continued Radial distribution of the relative excess of galaxies (n) with respect to the QSOs (in arcsec) for 
each quasar field and averaged over all the fields (indicated as mean). The thick and thin solid lines correspond to the 
objects detected in the i?-band and i?-band images, respectively. 



(1991) and Hall & Green (1998). Hintzen et al. (1991) 
found that their excess galaxies are ~0.5-0.8 mag brighter 
than first-ranked galaxies at this redshift, but only if no- 
evolution is assumed. This result has been used as an ar- 
gument for the presence of gravitational lensing in these 
objects (Webster 1991). However, as it has been quoted 
above, the observed galaxies could be at the redshifts of 
the QSOs if we take in account a small amount of evolu- 
tion, like in a E-type galaxy. 

Hall & Green (1998) found that the excess increases 
with power and spectral index of radio-emission. This ten- 
dency with spectral index was also recently reported by 
Mendes de Oliveira et al. (1998). Our results (Sect. 3.1) 



agree with the first tendency, although we can not test the 
second one due to the small baseline in radio spectral index 
of our sample ((a^g 6 )=-0.91±0.25). Hall & Green (1998) 
also found that the excess was distributed in two different 
components: a peak within 9 <40" from the quasars, and 
a smooth component along all the covered area. It is neces- 
sary to recall here that their results were based on images 
that allowed to sample the overdensities within 9 <^100" 
from the quasars. In Fig. 2, it is possible to see that in 
almost all our fields the overdensities present a weak peak 
within 9 <40" from the quasar, and a smooth component 
down to 9 ~160". This structure is present in the aver- 
aged distribution (Fig.2). Hall & Green(1998) found that 
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Fig. 3. Significance of the relative excess vs. radio power 
of the quasar. The dotted-line shows the result of a linear 
least-square fit over the mean significance. 

the peak is significant only in —25% of the fields, although 
the large-scale overdensity is present in all of them. This 
is also seen in our sample: only 1148+387 and 1206+439B 
fields present a strong peak in their central region. 

Three different scenarios could explain these results: 
(1) the quasars inhabit clusters that reside in large-scale 
overdensities of galaxies, as has been proposed by Hall & 
Green (1998), (2) they inhabit dynamically young, still 
not virialized, clusters, which expand smoothly to larger- 
scales (Ellingson et al. 1991 and Hall & Green 1998), and 
(3) not all the quasars are at the center of the observed 
overdensities. This last possibility has not been tested be- 
fore, and it could easily explain the observed distribution. 
In fact, only when the overdensity presents a clear peak in 
the inner region (within 6 <40") it is possible to say that 
the quasar is roughly at the center of the overdensity. 

In order to test this possibility we have estimated the 
density distribution of galaxies of our i?-band fields. The 
density distribution has been calculated using the gener- 
alized variable kernel estimate (Silverman 1986) using the 
Epanechnikov kernel (Epanechnikov 1969). This method is 
an improvement of the nearest neighbour approach, which 
adapts the amount of smoothing to the local density of 
galaxy counts. While the most widely used estimator of 
the density is based on the number of observations falling 
in a box of fixed width centred at the point of interest, 
the kernel estimator is inversely proportional to the size 
of the box needed to contain a given number of observa- 
tions (fc), and also inversely proportional to a smoothing 
parameter (h). The number of observations, k, were se- 



lected to be roughly the mean density of galaxy counts, 
and the smoothing parameter was arbitrarily fixed to two 
points. This density estimate defines better the properties 
of the distribution than the generally used (see discussion 
in Silverman 1986). 

Fig. 4 shows a contour plot of the density distribu- 
tion of galaxy counts for each field. The first contour was 
plotted at the level of the mean density, and the sepa- 
ration between sucessive contours is la. The position of 
the quasars in each field was marked as a blue star. It is 
clearly seen that the quasars are not at the maximum of 
the distributions, hence, they are not at the "center" of the 
detected overdensities. The distance between the quasars 
and the peak of the overdensity range between —40" and 
-100". 

We have recomputed the radial distributions of the 
relative excess taking as the origin the "center" of the ob- 
served galaxy densities. The result is shown in Fig. 5 for 
all the quasar fields and for the mean, where the position 
of the QSOs is indicated with an arrow. It is seen that 
these excess distributions are more intense at the inner 
region than the distributions around the quasar positions, 
and the overdensities within 35" are about twice more sig- 
nificant. In fact, there are fields (e.g. 1435+383) where the 
smoothness due to the off-centering of the quasar is dra- 
matic (compare Fig. 2 and Fig. 5). In all the fields (except 
perhaps 1435+383), the quasar is within the overdensity 
range. 

In Fig. 5 it is also seen that the excess still presents a 
smooth component in some of the quasars fields, although 
it is less significant than in Fig. 2. Therefore, although the 
off-centering artificially smooths the radial distribution of 
relative excess, our result does not completely exclude the 
other two hypothesis to explain the large scale or smooth 
component. 

5. Summary and conclusions 

In this study, the environments of 7 radio-loud quasars 
at a redshift range between ~1 and — 1.6 have been ob- 
served in the B and R band. A significant excess (up 
to — 3er) of faint galaxies (R >22, B >22.5) have been 
found around them, with magnitudes, angular size of 
the excess, and number of galaxies all compatible with 
clusters at the quasar redshifts. This result agrees with 
previous results on radio-loud quasars at similar red- 
shift ranges (Hintzen et al. 1991, Hall & Green 1998), 
and it is in contrast with results reported on radio-quiet 
ones (Boyle & Couch 1993). As at lower redshift ranges, 
the environments of radio-loud and radio-quiet sources 
appear to be signific antly different ( Yee fc Green 1987 , 
Ellingson et al. 199l[ ). 

We have found a tendency between radio power and 
excess of galaxies around QSOs. This trend could be ex- 
plained by two different mechanisms: (1) If the intracluster 
gas plays the role of the fuel of the nuclear engine, and ra- 
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Fig. 4. Contour plots of the estimated density distribution of galaxy counts in each field. The first contour corresponds 
to the mean density, and the separation between sucessive contours is la. The position of the quasar is marked with 
a blue star. The distance between small ticks is 20 arcsec. 



dio power is a good estimator of its activity, it is expected 
that radio power increases with gas density. It is also ex- 
pected that the amount of intracluster gas were larger 
in high density clusters, which could explain the correla- 
tion; (2) The steep-spectrum extended radio emission is 
explained by the interaction of jets with the intergalactic 
medium (e.g. Miley 1980). This medium, with a density 
of at least p > 10 -28 g- cm -3 , confines the radio lobes. 
It is expected that intergalactic medium were more dense 
in more populated clusters, and, therefore, there were a 
connection between cluster population and radio power. 
Maybe the combination of both mechanisms has to be in- 
voked to explain the proposed tendency. It is interesting to 



note here that at least the second hypothesis could explain 
the connection between excess and slope of radio emission 
found by Hall & Green (1998) and Mendes de Oliveira et 
al. (1998). 

Two possible explanations to the observed link be- 
tween nuclear activity and environment have been previ- 
ously proposed: the dynamical state of the cluster galaxies, 



the evolution of intracluster gas (Yee & Ellingson 1993). 
In both schemes the merger or interaction between the 
host galaxy and companions galaxies plays a substantial 
role, which reinforces the suggestion that the nuclear ac- 
tivity is triggered by merging processes (Barnes & Hern- 
quist 1992, Shlosman 1994). There are other evidences 
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Fig. 4. Continued Contour plots of the estimated density distribution of galaxy counts in each field. The first contour 
corresponds to the mean density, and the separation between sucessive contours is la. The position of the quasar is 
marked with a blue star. The distance between small ticks is 20 arcsec. 



that support this suggestion, like the large fraction of 
host-galaxies undergoing tidal interactions/merging pro- 
cess, their large absolute magnitudes, and the velocity 
distribution of galaxies around quasars (Smith et al. 1986; 
Hutchings 1987; Disney et al. 1995; Hutchings & Neff 1997; 
Carballo et al. 1998). Another possible explanation could 
be a connection between the mass of the central super- 
massive black hole and the depth of the potential well of 
the cluster, which could be reflected in the richness of the 
cluster itself. It is not clear which could be the primary 
physical cause among the above refereed, and maybe a 
combination of all of them would be needed to explain 
the results. 



We have found that the quasars of our sample are not 
located at the center of the clusters. The off-centering of 
the quasars partially explains the smoothing of the ra- 
dial distribution of the excess, which apparently shows 
a two-component distribution (Hall & Green 1998). A 
better determination of the nature of the overdensity is 
needed before to present a conclusive explanation for this 
off-centering. However, this result is not completly un- 
expected in the framework of merging as the trigger- 
ing/fuelling mechanism of the nuclear activity: Aarseth & 
Fall (1980) have shown that the merging process requires 
that the encounter velocity does not significantly exceed 
the internal velocity dispersion of the galaxies, approxi- 
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Fig. 5. Radial distribution of the relative excess of galaxies (n) with respect to the position of the center of the density 
(in arcsec) for each quasar field and averaged over all the fields (indicated as mean). The thick and thin solid lines 
correspond to the objects detected in the i?-band and i?-band images, respectively . The arrows mark the position of 
the quasar in the field. 



mately 200 km s . Therefore, in virialized cluster cores, 
galaxy-galaxy merging would be less effective that in the 
outer part of the clusters, where we have found that the 
quasars of our sample inhabit. 

Among the follow-up research projects suggested by 
this work, the most important is to select a sample of can- 
didates to cluster members among the galaxies detected 
in each field. It would be interesting to obtain optical-NIR 
multiband imaging and multislit spectroscopy of these 
candidates in order to verify the existence of overdensi- 
ties at the quasars redshifts, to discriminate the age and 
metallicity of the selected galaxies, and to do a dynamical 



study of the putative clusters in order to determine why 
the active nuclei are not at the center of them. 
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